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ABSTRACT 

The emission mechanism of TeV blazar 3C 66A is studied employing multi- wavelength 
observational data. The GeV spectrum observed by Fermi-LAT cannot be explained 
by pure synchrotron self-Compton (SSC) radiation. A multi-component leptonic jet 
model including the absorption of extragalactic background light reproduces the quasi- 
simultaneous spectral energy distribution. The external Compton component of the 
broad line region (BLR) is important for emission above 10 GeV, while emission from 
0.1 GeV to 10 GeV is dominated by the SSC component. Our model implies that 
the redshift of 3C 66A should be between 0.15 and 0.31, the most likely value being 
0.21. Our results also indicate that the 7-ray emission region is beyond the inner zone 
of BLR. The effects of model parameters and assumptions on the estimating of the 
redshift are discussed. 
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1 INTRODUCTION 

Blazars are the most extreme class of active galactic nu- 
clei (AGN). Their spectral energy distributions (SEDs) are 
characterized by two distinct bumps. The low-energy com- 
ponent originates in a relativistic jet, and consists of rela- 
tivistic electron synchrotron emission. The high-energy com- 
ponent could be produced by inverse Compto n (IC) scatter- 
ing (for a recent review, see lBottcherll2007h . Various soft 
photon sources seed the SSC and EC processes in the jet to 
produce gamma rays. 

Very high energy (VHE) photons emitted by extra- 
galactic sources are absorbed, through the pair-production 
process, by interaction with extragalactic background light 
(EBL) i|Stecker et al.lll992l ). The rate of absorption is de- 
termined by the EBL photon density and the redshift of 
the VHE source. The energy range of interest for back- 
ground photons is from optical to ultraviolet (UV). Since 
it is hard to measure the EBL directly, researchers have 
proposed three kinds of spectral model s: low limit (e.g., 
Kneiske et al.ll2010l: iRazzaque et al.ll2009 | ). mean level (e.g.. 
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Finke et al.ll201 0a; Franccschini ct al. 2008), and high level 

(e.g.. lStecker et al 1120061 ) . I Aharonian et al.1 (|2006h discussed 



some 7-ray blazars with unexpectedly hard spectra at large 
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redshifts, and sug gested that EBL is of the first type. 
lAlbert et ail l|2008n found th at the universe is more trans- 
parent t o 7-rays. However, IStecker et al.l l|2009l ) pointed 
out that lAlbert et al.l l|2008l ) do not significantly constrain 
the intergalactic low energy photon spectra and their high 
level EBL model is still valid. In an analysis of photons 
above 10 G e V from 7-ray sources detected by Fermi-LAT, 
IXbd o et al.l (|2010al ) found evidence to exclude the high 
level EBL model. By taking into account multi-frequency 
data ranging from radio to GeV, the SED of a blazar can 
be extrapolated into the VHE region. The redshift of the 
VHE source can then be constrained by comparing the de- 
absorbed VHE spectrum with the extrapolated one. The ac- 
curate determination of such redshifts helps us understand 
both the nature of EBL and the intrinsic emission model of 
VHE sources. 

3C 66A is classified as intermediate BL Lac (IBL), be- 
cau se its synchrotron peak lies between 10 15 Hz and 10 16 
Hz l|Perri et al.ll2003l : lAbdo et al.ll2010bl ). The most widely 
used redshift for 3 C 66A is 0.444, ba sed on a si ngle emission 
line m easurement l|Miller et al.lll978l ). However, iMiller et al.l 
(1 19781 ) stated that they were not sure of the reality of this 
emission featur e, and w a rned that the redshift was not 
reliable. Later, lLanzettal (|l993l ) confirmed the redshift of 
0.444 based on d ata from Inter n ationa l Ultraviolet Explorer 
(IUE). However, iBramel et ail (|2005l ) argued that the 3C 
66A redshift deter mined using IUE data is questionable. 
iFinke et all (|2008ah placed a lower limit on the redshift of 3C 
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66A, z ^ 0.096, using inf o rmati on regarding its host galaxy. 
Recently, IPrandini et al.l l|2010] ) concluded that the redshift 
of 3C 66A should be below 0.34 ± 0.05, and that the most 
likely redshift is 0.21 ± 0.05. These results assume that the 
de-absorbed V HE spectrum are n o t hard er than the Fermi- 
LAT spectrum. Ijoshi fc BottcheJ |2007l ) argued that 7-ray 
emission in the flare sta te could be dominated by an external 
Compton (EC) process. lYang fc Wanel i|2010h found that the 
VHE emission is contributed by EC when taking z = 0.444, 
or by pure SSC when z = 0.1. No te that this last res ult does 
not employ detailed Fermi data. lAbdo et all (|201ll ) fit well 
the SED of 3C 66A at flare state with the SSC+EC model, 
and suggested that the redshift of 3C 66A may be between 
0.2 and 0.3. 

A quasi-simultaneous, multi-wavelength campaign of 
observations for 3C 66A was carried out by Fermi and Swift 
from August 2008 to November 2008. Data from the ra- 
dio, optical, U V, X-ray, and GeV 7-ray to TeV 7-ray bands 
are available i|Abdo et all l2010bl ). VERITAS observed 3C 
66A fr om September through November 2008 (|Acciari et al.l 
1200931 ) . Assuming the redshift of 3C 66A is 0.444, its VHE 
spectrum corrected using the mean le vel EBL model is 
harder than the obse r ved Ge V spectrum l|Abdo et al.l l2010b; 
lAcciari et all l2009al . 1201031 ). In this work, we fit the up- 
to-date SED of 3C 66A with a multi-component lepton jet 
model that takes into account EBL absorption. This model 
constrains the redshift to lie between 0.15 and 0.31. In Sec- 
tion 2, we describe the model and constrain its free param- 
eters. Our results are presented in Section 3. A discussion 
of the data and our conclusions are given in Section 4. We 
adopt the cosmological parameters (Ho, £im, ^a) = (70 km 
s" 1 Mpc" 1 , 0.3, 0.7) throughout. 



2 MODEL AND PARAMETERS' 
CONSTRAINING 

We assume that 7-ray emission is produced in a spheri- 
cal blob of the jet, which is moving relativistically at a 
small angle to our line of sight. The observed radiation is 
strongly boosted by a relativistic Doppler factor So- The 
relativistic electrons inside the blob lose energy via syn- 
chrotro n emission and IC s cattering. The electron distribu- 
tion is (iDermer et al.ll2009l) . 



^0(7') = K' e H(-y'; 7^,7^)7' P1 exp(--y' / '7b) 

xH[( P2 - Pl ) 7 ; - 7'] + [(pa - Pl Hr- pi -y'- p2 

0'min 1 7max 

), (l) 



where H(x;xi,X2) is the Heaviside function: H(x;xi,X2) = 
1 for xi ^ x ^ X2 and H(x; xi, X2) = everywhere else; 
as well as H(x) — for x < and H(x) = 1 for x ^ 0. 
In the co-moving frame, this distribution is a double power 
law with two energy cutoffs: 7mi n and 7m ax . The spectrum is 
smoothly connected, with indices p\ and p2 below and above 
the electrons' break energy 7 b . Note that here and through- 
out the paper, unprimed quantities refer to the observer's 

frame and primed ones refer to the co -moving frame . 

The multi-component model of IDermer et all (|2009l ) 
is used to calculate the multi-wavelength SED. This 
model takes into account 7-ray attenuation by the BLR- 
scattered radiation field. For EC components, we con- 



sider photons emitted directly from the accretion disk, us- 
ing BLR-scattered photons as the seed photons. Here we 
brie fly describe th e mod el for numerical calculation, please 
sec 1 IDermer et ail (|2009l ) in detail. The vF v spectrum of 
Compton-scattered external photons is 
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where 7 = 7'fo, d^(z) is the luminosity distance, z is the 
source redshift, c is the speed of light, fj,* is the collision an- 
gle between the relativistic electrons and seed photons and 
w(e*,fi*) is energy density of external photons with energy 
e». The Compton cross section in the head-on approxima- 
tion, accurate for relativistic electrons from the Thomson 
through the Klein-Nishina regime, is 
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where r c is the classical electron radius, the invariant colli- 
sion energy e = 7e*(l — fi,) and y = 1 — e s /7- The lower 
limit on the electron energy 7i ow and the upper limit e*, up 
are 
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Gamma-ray photons are subject to 77 attenuation by 
the ambient photons in the environment of the black hole. 
The absorption probability per unit pathlength is 

i^_internal { 



dl 



dfi*(l - ft,) 



w(e», f2«) 



,{P)de, 



(6) 

l|Dermer et al.ll2009l ). where fj, t is the angle between the 
direction of the Gamma-ray photon with energy ei and 
the target photon with energy e, and m e is the rest mass 
of electron(the scattered photon travels in the same direc- 
tion as t he relativistic elect r on). T he cross-section a- ry (0) is 
given bv lGould fc Shrederl l| 19671 ), with f3 = ^/l - Ifs and 
s = eie*(l + z)(l-jt,)/2. 

We assume that the BLR is a spherically symmetric 
shell with inner radius Ri and outer radius R . We also as- 
sume that the gas density of the BLR has the power-law 
distribution n c (r) = n>o(jr) > where Ri ^ r ^ R . The ra- 
dial Thomson depth is given by tt = ctt Jn ° drn e (r), where 
r is t he distance from the central black hole (|Dermer et alj 
2009). In our calculation we use tt = 0.01, wh ich is 
the typical value (e.g., iFinke et ail l2010bl: iReimerl 120071 : 
iDonea fc Prothero j2003 ). lKaspi fc Netzed (|l999l ) suggested 
that the particle density of BLR scales as r _1 '° or r -1,5 . In 
our calculation, we adopt the expo nent C = — 1-0 . 

Using reverberation mapping, iBentz et alj (|2009h de- 
rived an improved empirical relationship between BLR ra- 
dius -Rblr and luminosity L\ at 5100A: 

log(i?BL R .) = -21.3 + 0.519 • log(AL A (5100A)) . (7) 



The V-band magnitude of 3C 66A 



15.21 



Emission and redshift of 3C 66A 3 



jVeron-Cettv fc Veron Il2010l). W e use the optical spectral 
index given by Fiorucci et al.l (|2004l ) to calculate the 
average flux at 5100A, which is 2.785 mj. In this work, we 
take the estimated .Rb lr as the outer radius of the BLR R . 
iPeterson etail l|l994 ) suggested that the typical size of the 
BLR in quasars is on t he order of lig h t-months. We follow 
sever al authors (e.g., iReime 3 120071 ; iDonea fc Protherod 
2003), using the relationship Ri = R o /A0 to derive a value 
for Ri. 

To simplify calculation, the BLR-scattered photon field 
is assumed to be monochromatic with ene rgy e«, which 
is th e mean energy from the accretion disk (|Dermer et al.l 
2009). The approximate mean dimensionless photon energy 
at radius R from a standard accretion disk is given by 
l|Finke et al.ll2010bl ) 



e d (R) = 1.5 x 10" 4 ( 



104 



R, 



M 8 T) 



(8) 



The accretion luminosity is ^Edd = r d i which here has 
the value 0.03. The Eddington luminosity is I/Edd = 1-26 x 
10 46 Mgergs ■ s~\ and L d is the accretion disk luminosity. 
The accretion efficiency r\ is 0.1. The gravitational radius 
is r g = Of 1.5 x 10 13 M 8 cm. T he mass of the 3C 66A 
black hole is M s = yff™^- = 4.0 (|Ghisellini et alj l2010h . 

The typical energy of pho tons from the disk is ed(2.25Rf n ) 
l)Shakura fc Sunvaevlll973l ). where Rf n is the inner radius of 
accretion disk. Rf n is approximately equal to 3r g for a Kerr 
black hole, or 6r s for a Schwarzschild black hole. In this 



work, we adopt e* = £d(10r g 



x 10 J . The energy 



density of BLR-scattered photon field is 

L d r 2 



i(e*,^i*;r b ) = 



3cr b 



(9) 



ijDermer et al.l l2009b where r b is the distance from the emis- 
sion blob to t he central black hole and F(fi t , r b ) is the func- 
tion given bv lDermer et all l|2009h (their Eq.(97)). 
The intrinsic flux from extragalactic sources is 
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where /observed is the measured VHE flux, and t(E 1 ,, 
the optical depth of a 7-ray with energy E~ ( at redshift z. 
T(E~f, z) is given by (|Franceschini et al]|2008 
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jj^ gives the relationship between cosmic 
time and redshift, which is = H (1 1 +Z1 ) [(1 + zi) 2 (^ m zi + 

l)-2i(zi+2)n A ]~ 1/2 . dn( ^ zi) is the number density of EBL 
photons with energy e at redshit Z\, and x = 1 — cosO, where 

9 is the collision angle. 

Here, we use the EBL model of iFranceschini et all 

i 2IKIs£] to de-absorb the observed VHE SED. This model 
is based on observations and takes into account all available 



1 Opacities for photon-photon interaction as a function 
of the source redshift are available on the the website 
http://www.astro.unipd.it/background 



Table 1. Model Parameters for Fig. [TJ 



parameters 


z = 0.15 


z = 0.21 


2 = 0.31 


7 b (10 3 ) 


5.8 


6.3 


7.6 


7 nlin (10 3 ) 


1.93 


1.90 


1.76 


&D 


38 


36 


43 




0.69 


1.17 


1.21 


Ro (pc) 


0.1 


0.14 


0.22 


r b (Ro) 


1.03 


0.89 


0.72 


B/Bcq 


0.306 


0.335 


0.352 


x 2 


3.099 


1.498 


3.880 



infor mation on cosmic sources c ontributing background pho- 
tons. IFranceschini et all l|2008l ) also considered the redshift 
evolution of EBL. 

Several p arameters in our mod el can be constrained by 
observations. iBottcher et all d2005 | ) estimated t h e Dop pler 
factor of 3C 66A as <5 D » 15. IBottcher et al.l (|2009h ex- 
cluded extreme values of the Doppler factor in the range 
So 50. The size of the emission blob can be constrained 
by the observed variability timescales t vaT , because i? b = 
tv,mm<5zjc/(l + 2) < 5j3ct t ar/(l + z). Here Rb is the radius 
of the blob in the comovi ng frame, and t v ,min is the small- 
est variability timescale. iTakalol (| 19961 ") reported a micro- 
variability with £ var ~ 2.16 x 10 4 s and Amag ~ 0.2. 
lAbdo et al.l <|201ll > reported shorter variability at optical 
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1.44 x 10*s. 



3 THE RESULTS 

From the observed photon index in the X-ray band Tx 
l|Abdo et al.ll2010bl ), we can roughly estimate that p2 ~ 4.0 
using the relationship Px — 1 = \ (p2 — 1) . In our calculation, 
the values of pi and p2 are fixed at 2.0 and 4.0 respectively. 
There are few constraints on the values of the magnetic field 
B and 7max. We use the fixed parameters B — 0.168 G and 
7max = 3.0 x 10 6 . Table[T]lists all other model parameters, in- 
cluding £?cq (the equipartition magnetic field). The reduced 
X 2 is calculated from the normalized residuals. 

In Fig. [T] we show the modeling results at different 
plausible redshifts. The filled circles are quasi-simultaneous 
data in bands ranging from radio to GeV. The observed 
VERITAS data are de-a bsorbed using the EBL model of 
IFranceschini et all (|2008l ) with different redshifts. We obtain 
the best fit when z = 0.21. It is found that the accretion- 
disk component is negligible compared to the SSC and BLR 
components. SSC and EC are responsible for emission in the 
GeV-TeV range. Emission from 0.1 GeV to 10 GeV is dom- 
inated by SSC. Above 10 GeV, the EC component of BLR 
is more important. 

In Fig. [2] we show the 77 absorption by BLR-scattered 
radiation at different locations of the blob when taking 
2 = 0.21. There is significant absorption when the blob is 
inside the inner zone of the BLR. Beyond the inner zone, 
absorption is negligible. 
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Figure 1. In panel (a), we show the best modeling results with 
z = 0.21. The filled square are the de-absorbed VHE data with 
z = 0.21. The dashed, dash-doted, dotted and thick solid lines are 
SSC component, accretion-disk, BLR-reproduced component and 
the sum of multi-component, respectively. In panel (b), the open 
square, filled square and open circle are the de-absorbed VHE 
data with z = 0.15, 0.21 and 0.31, respectively. The dash-dotted, 
solid and dashed lines are the model results at z = 0.15, 0.21 an d 
0.31, respectively. All observed data are from lAbdo et al] fcOlObl) . 
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Figure 2. 77 optical depth for 7-ray interaction with BLR- 
reproduced photons at different distances from central BH when 
z = 0.21. 



4 DISCUSSION AND CONCLUSIONS 

At present, three IBLs (3C 66A, W com and PKS 1424+240) 
have been detected in VHE region. The first two require 
an EC component to model high-energy emiss ion (e.g., 
al] l2009bl . l2010bl ; lYang fc Wanel |2010h . A pure 
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SSC model fails to explain the GeV spectrum of 3C 66A 
observed by Fermi-LAT. However, a satisfactory fit to these 
data can be obtained by our multi-component model (see 

Fig. mi. 

The multi-component mo del requires a la r ge va lue of 
7^ in ~ 2 x 10 3 . As argued by iTavecchio et al.l jjooj ), this 
outcome seems to provide important clues to the electron 



acceleration process and the role of energy loss. A large value 
of j' miD leads to a steep spectrum in the low-energy band, so 
our model does not explain the observed radio emission. This 
may come from a larger emission region. Assuming the total 
particle power is te n times of the electron power in the jet 
l|Finke et alj|2008bh . the ratio between the best-fit magnetic 
field and equipartition magnetic field (see Table [TJ is close 
to a state of equipartition between the total particle energy 
density and the magnetic field energy density. 

A simple assumption is that the BLR scattered pho- 
ton field is characterize d by an isotropic bla ckbody (e.g., 
Ilnoue fc Takaharal [l99r3 ; lYang fc Wanel I2010T ). Our model 
takes into account not only the specific shell structure of the 
BLR, but also the angular dependence of the photon distri- 
bution. As shown in panel (a) of Fig. [1] BLR-sca t tered pho- 
tons are important to the EC process, fran et all (|2006h also 
found that soft seed photons for EC predominantly come 
from the BLR in a sample of 7-ray loud blazars. Further- 
more, our results show that the contributions of different EC 
components to the 7-ray spectrum are dependent on the lo- 
cation of the blob. As the distance between the central BH 
and the blob increases, the BLR component grows more im- 
portant relative to the accretion-disk component. The BLR- 
scattered radiation field not only provides the seed photons 
for IC, but also attenuates the 7-rays. In flat-spectrum radio 
quasars (FSRQs), absorption by BLR radiatio n is significant 
and can cause a break in t he GeV band (e . g.. lReimerll2007l; 



Sitar ek fc Bednarekl 120081 ; iBai et all 120091 ; iPoutanen et al.l 
Our results show that 7-ray absorption by BLR ra- 



diation becomes more and more slight as the value of rb 
increases (see Fig. [5J. 

Thanks to these high-quality, multi-frequency data, we 
can use the extrapolated intrinsic VHE spectra to constrain 
the redshift of 3C 66A. As shown in panel (b) of Fig. [T] the 
de-absorbed VHE spectrum is too steep to be fitted well by 
our model if the redshift is below 0.15. However, the same 
is true if the redshift is above 0.31, as the shape of the de- 
absorbed VHE spectra becomes strange. Hence, the redshift 
of 3C 66A should be between 0.15 and 0.31. We calculate 
the reduced % 2 f° r H 7 _ra y spectra at different redshifts 
between 0.15 and 0.31. The minimum value, x 2 = 1-498, oc- 
curs for z = 0. 21, which mean s the m ost likely redshift for 
3C 66 A is 0.21. |Prandini et all (|2010h assume that the VHE 
spectrum corrected for the absorption through the interac- 
tion with the EBL can not be harder than the spectrum in 
the Fermi/LAT band. They derived a EBL-corrected spec- 
tra as a function of the redshift with a TeV blazars sample of 
known redshift. Comparing the TeV corrected spectra with 
GeV LAT spectra, they get the estimation distance of un- 
known redshift blazars. They suggested that the most likely 
redshift of 3C 66A is 0.21 ± 0.05. Two results from different 
methods are completely consistent. 

The contribution of the BLR component is crucial for 
constraining the redshift of 3C 66 A. The BLR structure (i?i, 
^o, Cj t t) and the characteristics of the central source (the 
black hole and its accretion disk) can affect the contribution 
of the BLR component. R can be constrained by Eq.(7). 
We assume (^ Edd , R /Ri, Tr, C) = (0.03, 40, 0.01, -1) and get 
the best model fitting with z — 0.21. The effects of these pa- 
rameters and assumptions on the estimating of the redshift 
are discussed. Fig. [3ja), (b), (c) and (d) show the effects of 
accretion luminosity, BLR size, radial Thomson depth and 
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Table 2. Model Parameters for Fig. [3] The blank ones represent 
these parameter are the same as that we adopted when z = 0.21. 



parameters 


fedd = 0.01 


f =5 


T T = 0.1 


C = -2 


| 


7 b (10 3 ) 

7^m(10 3 ) 

So 

t 1 ,,min(10 4 s) 


1.8 
1.2 


2.0 


5.6 
2.5 
37 
1.05 


3.2 




r h {Ro) 


0.65 


1.02 


1.31 


0.52 




B/Beq 


0.333 


0.344 


0.366 


0.794 




x 2 


1.251 


1.511 


1.934 


1.314 


E 



accretion luminosity=0.01 




distribution exponent of BLR gas density on the model fit- 
ting respectively. For clarity, only the modeling results in 
the high-frequency range of the case of z = 0.21 are shown. 
The model parameters are listed in TableEJ The parameters 
which we don't particularly claim are the same as that we 
used in the case of z = 0.21. The 77 absorption by BLR- 
scattered radiation is also negligible in all of the four cases. 

From the reduced \ 2 ( see Table [2]), we find the fitting 
results are in agreement with that we originally obtained. 
We therefore conclude that the assumptions of BLR struc- 
ture and the characteristics of the central source don't af- 
fect our estimating on the redshift. It can be interpreted as 
that these assumptions affect the density of external pho- 
tons field, and we can adjust the location of the blob (Yb) 
to determine the suitable density of external photons field 
to reproduce the GeV radiation. Therefore, determining the 
location of the blob is important for constraining the struc- 
ture of the BLR. It's found that the reduced \ 2 gets better 
when iEdd = 0.01 and worse when tt = 0.1 compared with 
the case of ^Edd = 0.03 and tt = 0.01, respectively. Hence, 
the more appropriate parameters may be ^Edd = 0.01 and 
tt = 0.01 for 3C 66A. From Table[2] we see that the changes 
of parameters are slight and reasonable except the case when 
( = —2 (see Table. [2J. In the case of ( = —2, the den- 
sity of the BLR-scattered photon field declines more quickly 
with the blob's distance from the central black hole. For 
fitting the GeV data, the smaller value of is required, 
then the accretion-disk component becomes more impor- 
tant. As shown in panel (d) in Fig. [3] radiation from 0.1- 
10 GeV is dominated by accretion-disk component. In this 
case, lager value of t v ,min = 3.2 x 10 4 s is required, which is 
in contrast with the observed optical variability time scale: 
ivar ~ 1-44 x 10 4 s. Our results therefore indicate that the 
case of C = —2 seems to be excluded for 3C 66A. We also 
conclude from model fitting that the blob is located beyond 
the inner zone of the BLR, may be between 0.5R o and 1.37? - 

In conclusion, we study the multi-frequency emission of 
3C 66A by using multi-wavelength observational data over 
radio through GeV-TeV energy region. Under the mean- 
limit EBL absorption, the model constrains the redshift of 
3C 66A to lie between 0.15 and 0.31. The most likely value 
is 0.21. 



Figure 3. The effects of different assumptions of BLR structure 
and the characteristics of central source on the estimating of the 
redshift. The symbols are the same as that in Fig. [T1 
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